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Introduction 

This is a final report on Grant No, NSG 4018. The project began on 
January 1, 1979 and was to have terminated on December 31, 1979. A no- 
cost time extension was requested and approved extending the termination 
date to June 30, 1980. 

Dr. Robert L. Swaim, Principal Investigator, has devoted ten percent 
time throughout the project. One Ph.D. student, Mr, Supat Poopaka, has 
been employed fifty percent time throughout the project. Mr. Poopaka is 
completing his Ph.D. thesis on the project and expects to complete it and 
take the final examination, where he defends the thesis, in September 1980. 
Three copies of the thesis will be submitted to the Technical Officer, Mr. 
Glenn Gilyard, and two copies to the NASA Scientific and Technical Infor- 
mation Facility in late September. 

As the details of the research will be documented in the thesis, 
th*s final report on the project will be just a brief summary. 

Problem Statement 

The effects on handling qualities of low frequency symmetric elastic 
mode interaction with the rigid body dynamics of a large flexible air- 
craft was analyzed by use of a mathematical pilot modeling computer 
simulation. An extension of the optimal control model for a human pilot 
was made so that the mode interaction effects on the pilot's control task 
could be assessed. Pilot ratings were determined for a longitudinal 
tracking task with parametric variations in the undamped natural frequencies 
of the two lowest frequency symmetric elast'r modes made to induce varying 
amounts of mode interaction. 

A summary of our mathematical formulation of the problem is contained 
in Attachment 1, which is a paper presented by Mr. Poopaka at the Eleventh 
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Southwestern Graduate Research Conference in Applied Mechanics held April 
11-12, 1980 at Oklahoma State University. 


Results 

Our approach of relating numerical performance index values associated 
with the frequency variations used in several dynamic cases to a numerical 
Cooper-Harper pilot rating has proved successful in discriminating when 
the mathematical pilot can or cannot separate rigid from elastic response 
in the tracking task. The detailed presentation of these results will be 
contained in the forthcoming Ph.D. thesis by Mr. Poopaka. An outline of 
the thesis content is presented below. 


Thesis Outll ne 


I. Introduction 

A. Handling Qualities and Pilot Ratings 

B. Pilot Rating Assessment Techniques 

C. Background on the Effects of Dynamic Aeroelasticity 

on Handling Qualities 

D. Objectives and Scope of Study 

E. Plan of Presentation 

II. Literature Survey 

A. Past Results on Elastic Airplane Research 

III. Equations of Motion 

A. Longitudinal Equations of Motion for a Flexible Airplane 

B. Turbulence Model 

C. Attitude Director Equations 

D. B-l Flight Condition 

IV. Pilot Modeling 

A. Manual Control and Pilot Models 

B. Optimal Control Model 

C. Internal Model 

D. Closed- Loop Performance Equations 

E. Hess's Pilot Rating Method 


Thesis Outline (cont.) 

V. Effects of Elastic Mode Interaction on Handling Qualities 

A. Define Controllable Boundary 

B. Controllable Boundary Based on Open-Loop Parameters 

VI. Conclusions and Recommendations 

VII. Selected Bibliography 

VIII. Appendix A. Numerical Values of Stability Derivatives 
and Equations of Motion 

IX. Appendix B. Derivations Related to Singular Perturbation 
Theory 

X. Appendix C. Computational Algorithms 


Project Reports and Presentations 

A list of the reports and presentations generated on the project 
is given below. 


1. R. L. Swaim and Wen-Yo Yen, "Effects of Dynamic Aeroelasticity on 
Aircraft Handling Qualities," Journal of Aircraft , Vol. 16, No. 9, 
Sept. 1979, pp. 635-637. Also presented at AIAA Oklahoma Section 
Mini -Symposium, University of Oklahoma, Feb. 24, 1979. 

2. R. L. Swaim, "Handling Qualities of Large Flexible Aircraft," 
presented at AIAA Oklahoma Section Mini -Symposium, Oklahoma State 
University, Feb. 23, 1980. 

3. Supat Poopaka, "Handling Qualities of Large Flexible Aircraft," 

Ph.D. Thesis Research Proposal, Oklahoma State University, Nov. 21, 
1979. 

4. Supat Poopaka, "Handling Qualities of Large Flexible Aircraft," 
presented at Eleventh Southwestern Graduate Research Conference in 
Applied Mechanics, Oklahoma State University, April 11-12, 1980. 

5. R. L. Swaim, "Handling Qualities of Large Flexible Control -Configured 
Aircraft," Semi-Annual Status Report on Grant No. NSG 4018, June 22, 
1979. 

6. R. L. Swaim, "Handling Qualities of Large Flexible Control -Configured 
Aircraft," Semi-Annual Status Report on Grant No. NSG 4018, Dec. 31, 
1979. 
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7. R. L. Swaim, "Handling Qualities of Large Flexible Control -Configured 
Aircraft," Final Technical Report on Grant No. NSG 4018, June 30, 1980. 

8. Supat Poopaka, "Handling Qualities of Large Flexible Aircraft," Ph.D. 
Thesis, Oklahoma State University, School of Mechanical and Aerospace 
Engineering, to appear Sept. 1980. 

9. Supat Poopaka and R. L. Swaim, "A Mathematical Pilot Model for Appli- 
cation to Highly Elastic Aircraft Handling Qualities Studies," under 
preparation for submission to the Journal of Guidance and Control in 
Sept. 1980. 
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The effects on handling qualities of elastic mode interaction with 
the rigid-body dynamics of a large flexible aircraft are discussed. An 
extension of the optimal control model for the human operator is made 
so that the mode interaction effects on-the pilot's control task can be 
assessed. 

Introduction 

The handling qualities are the flying qualities of a piloted air- 
craft. They are the characteristics of an aircraft that govern the 
ease and precision with which a pilot is able to perform the control task 
required in support of the aircraft mission flight phase. The pilot's 
subjective opinion on the handling qualities is called a pilot rating 
and usually based on the Cooper-Hjper pilot rating scale 11 J. 

The effects of elastic mode interaction on the handling qualities 
have been reported recently in Reference 2 based on a ground-based pilot- 
flown simulation. In this report we develop an extension of the optimal 
control model for the human operator so that the mode interaction effects 
on handling qualities can be easily assessed. 

Problem Statement 

The dynamics of an aircraft to be controlled by the human pilot are 
described by a set of small perturbation equations of motion: 

'K a S t '> • + 6*.a<v; f 

y«(i) - Ca.Xcd) + DU.fi) 


( 1 ) 
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where x, 8 col. lx,„,x ai ,,x_ J , » , Is a turbulence state vector, x a „ 

d dy ar uU itij ai 

is a rigid body states, x 3 _ is an elastic body states, u is a control 
input, y. is a displayed output, and w is a zero-mean, gaussian, white 

O 

noise process with autocovariance. 

E{ vCt,) * W <U± ( -t 2 ) 

For a cruise, level flight condition, eqn. (1) represents the 
longitudinal equations of motion and y '(t) is an attitude director equa" 

a 

tion. 

1L 

where x p indicates pilot fuselage station, <f>j(x p ) the slop^ of the j n 
symmetric elastic mode at that station, £-(t) the generalized displace- 

* J 

ment,' 0(t) the rigid body pitch angle, and Q a (t) the elastic contribution 

6 

to total pilch angle (see Fig. 1). 



The task of the pilot is to keep the rigid pitch response following 
the command director. Since only the total pitch response is shown to 
the pilot, therefore he has to distinguish rigid body motion from the 
total motion. We want to establish the boundary between when the pilot 
can visually separate the rigid body /notion from|;he total motion in terms 
of undamped natural frequencies of the elastic modes. The rigid body 
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mode parameters will he maintained at the values knuwn to (jive good 
handling qualities by* the pole- placement technique while the elastic 
mode frequencies will be parametrically varied to introduce the mode in- 
teraction. The pilot model developed for computer simulation of this 
task is presented in the next section. 

The Pilot Model 

The optimal control model for the human operator developed here is 
a modified version of the standard optimal control model (3,4). The model 
structure is illustrated in Fig. 2. 



Fig. 2 Optimal Control Model for Human Operator 

The relation u(t) = r(t-x) is approximated by a first-order Padd’* 
approximation which can be expressed in the state variable form as 

ulc-o - •zcv) - rcO 

icy) - r ary; + | r **r/) ( 4 ) 

The perceived information y(t) is a noisy version of y (t), i.e., 
- P a 


( 5 ) 
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where v y (t) Is a zero-moan, gau^.i.in, white noise with autocovariance 
E{vya t )iryrti)j * V, ( 6 ) 

and V u is known to scale with mean square of y a ; 
y « 

V B . « **»{*/«*] (7) 

where the observation noise/signal ratio e B O.Oln. 

y i 

The pilot generates a command control input m(t) which is then 
transformed to r(t) via the relation 

rco + ra) ~ wa) ^ 


where T n is the neuromotor lag matrix and v n) (t) is a zero-mean, gaussian, 
white noise, with autocovariance 

£ £ j <5 Urt t ) 

2 

and V m is known to scale with E{m j ( t ) ), i.e., 

V*. = ' ( 10 ) 


where the motor noise/signal ratio e m = 0.003n. 

Equations (4) and (8) may be augmented to (1) to define an aug- 
mented system of equations 


% c CV = A C X C U) 4 & c mCV +E c rt c CO 
k) a (+) « C c 'XcCV 
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The aircraft dynamics of eqn, (1) can be written in the form 
%'(+) ■** fin K\l*\ "f A'*- B/ <A(i) ■+£,W(i') 

* Aafi.a) 4- fiit -f Bi uii ) +£»^ 

( 12 ) 

* £/ (i) 4- C 2 Ki.fO + D U.C4) 
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where x.| is a vector of the rigid body state variables, Xg is a vector 
of the elastic mode state variables, and \t is a small positive scalar 
constant which can be unknown in this study. 

In performing estimation and control command generation, it is nec- 
essary for the pilot to have a knowledge of the aircraft dynamics, the 
human limitation parameters {t ,T n »Vy , V m > and weighting coefficients of the 
performance index. This knowledge is called the internal model. In our 
aircraft dynamics model, the display consists of a slowly varying part 
due mainly to the rigid body dynamics and a high frequency oscillation 
part from the elastic mode. From the past experiment (2 J it is evident 
that the pilot ignores the low amplitude high freuqency oscillation part 
of the display. Therefore, in the pilot model, we will use the slowly 
varying dynamics subsystem as the internal model. 

By the singular perturbation technique [5,6], we can decompose the 
slowly varying part from the system (12) by letting p-K> + , from which we 
get 

= fid. 

y d (4) * u«> (13) 
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The pilot control task is assumed to be adequately reflected in the 
choice of a control r(«) that minimizes the performance index 

JOO-jJ V <; + (14) 

*■ 

conditioned on the perceived information y (•)• The command control m(t) 

r 

of eqn. (8) is then given by 

v*(4)e~L or i$j(i) (15) 


where 
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The weighting Q r is chosen such that the resulting T n is equal to 0.1 
second, which is a typical value for the human neuromotor lag. 

The state x t ( t ) is the best estimate of x t (t) generated by a Kalman 

t 

filter, 

£ t Cf)= A±X*Wf B 4 *»«>■*■ (17) 
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where D satisfies the equation 


O s A t Z + lo W L 



( 18 ) 


Combining eqns. (5) with (11), (15) and (17) yields the closed loop 
system 

AcVcd) ^ € C L* +E c * c H> 

% t «) - + 2; c," V a " ' [ Q - C. x*. * ly oj (19) 


The covariances E(y a (t)**) and E{m.{t)^l for eqns. (7) and (10) can be 
obtained by solving eqns. (19) for cov. (x c ,x t ). 

Conclusion 

The technique developed here is being applied to determine the 
handling qualities of a large flexible aircraft of Reference 2 in a longi- 
tudinal tracking task with parametric variations in the undamped natural 
frequencies of the two lowest frequency, symmetric elastic modes made to 
induce varying amounts of mode interaction. The pilot rating is found by 
using the technique of Reference 7. 
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